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Abstract 

Theory of exciton fine structure in semiconductor quantum dots and its dependence on quantum 
dot anisotropy and external lateral electric field is presented. The effective exciton Hamiltonian 
including long range electron-hole exchange interaction is derived within the k * p effective mass 
approximation (EMA). The exchange matrix elements of the Hamiltonian are expressed explicitly 
in terms of electron and hole envelope functions. The matrix element responsible for the "bright" 
exciton splitting is identified and analyzed. An excitonic fine structure for a model quantum 
dot with quasi- two-dimensional anisotropic harmonic oscillator (2DLAHO) confining potential is 
analyzed as a function of the shape anisotropy, size and applied lateral electric field. 
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I. INTRODUCTION 



[l[ [2} of semiconductor quantum dots (QDs) 0, Q] 



One of the promising applications [l|, |2j of semiconductor quantum dots (QDs) [3|, |4| for 
uantum cryptography is the generation of entangled photon pairs (EPPs) on demand 



a a. 



15. 



EPPs can be generated with great efficiency via the biexciton cascade process 
(BCP) [l| in which a biexciton radiatively decays into the ground state via two different 
indistinguishable paths involving two intermediate dipole-active ("bright") exciton states. 
The major barrier to EPPs generation in BCP is the splitting of the intermediate "bright" 
exciton levels which distinguishes the paths of radiative decay and, as a result, destroys the 
entanglement. The splitting and mixing of the two bright exciton states is controlled by 
the long ranged electron-hole exchange (LRE) interaction and depends on dot asymmetry 
and applied magnetic field 0, Lateral electric fields have been applied in the hope to 



ff 



12|. 



control the dot anisotropy and hence anisotropic exchange splitting [lOj, 

Better understanding of the electron-hole exchange interaction in QDs, particularly its 
LRE part, should help in the development of EPP generation schemes. 

The exchange integral can be decomposed in real space into long-range part, i.e., exchange 
interaction between two "transition densities" localized in two different Wigner-Seitz (WS) 
cells and short-range part, i.e. exchange interaction within a single WS cell. A closely related 
decomposition of exchange into analytical and non- analytical part exists in the reciprocal 
space. In what follows, we will use the "real space" definition of long-range and short-range 
exchange (SRE). We will reserve the use of words "local exchange" and "nonlocal exchange" 
to describe the type of integrals which contribute to the exchange matrix elements. 

LRE electron-hole interaction in bulk semiconductors was investigated almost 40 years 
ago 
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15] and it is now well established that LRE splits the energy levels of "bright" 



excitons. For example, in zinc-blende direct band gap semiconductors with four-fold de- 
generate valence band r 8v and two-fold degenerate conduction band r 6c , the ground state 
exciton is eight-fold degenerate. The addition of SRE interaction into the excitonic Hamil- 
tonian will split the eight-fold degenerate ground state into "dark" and "bright" multiplets 
with degeneracies of five and three. The addition of LRE interaction will further modify the 
fine structure by splitting three-fold degenerate "bright" exciton level into two transverse 
excitons with J z = ±1 and one longitudinal exciton with J z = . 

Recent advances in single QD spectroscopies motivated re-examination of electron-hole 
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exchange in systems with reduced dimensionality within the framework of envelope function 



approximation 



tfi, 12, 3 



19 



20 



21 



22j and, more recently, from the point of view of 
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26j. 



161 ] have shown that if the 



empirical atomistic tight-binding approximation 

Within the envelope function approximation, Takagahara 
electron-hole pair envelope contains only Y 00 (6, 0) angular momentum component, then the 
long-range part of the electron-hole exchange vanishes. In Refs. \v\ , [lil . fine structure of 
localized exciton levels in quantum wells was considered. Efros el. al. 18| and Gupalov et. 
investigated band-edge excitonic fine structure of spherical CdSe nanocrystals. In 
20H . Takahagara derived effective eight-band excitonic Hamiltonian which takes into 




account electron- 
malism of Refs. 



no 
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e exchange interaction. In particular, within the envelo pe function for- 



201 ]. LRE is a dipole-dipole interaction. Takagahara [20j applied his 



scheme to investigate excitonic fine structure of disk-like GaAs/AlGaAs QDs as a function 
of QD anisotropy and size. It was demonstrated numerically in Ref. [20] that LRE vanishes 
for the ground state "bright" exciton doublet in QDs with rotationally symmetric 3D confin- 
ing potential. A similar to that of Ref. 20j representation of the electron-hole exchange was 
discussed by Maialle et. al. [lTj] in connection with the exciton spin dynamics in quantum 
wells. 

Within atomistic empirical tight-binding approach it has been demonstrated that the 
physical origin of long-range exchange interaction might be different from that of bulk 3 
and that LRE has a nonvanishing magnitude even in "shape-symmetric" dots 25|, |26 |. 
The former and the latter stem from the loss of local orthogonality on the unit cell scale 
between electron and hole single-particle orbitals compared to the bulk case and "reduced" 
symmetry of the atomistic confining potential, respectively. For T 6 x T 7 exciton, Gupalov 



et. al. 



24J identified the dipole-dipole and monopole-monopole contributions to the LRE 



with the intra-atomic and interatomic transition densities, respectively. 

In this work, we will re-examine electron-hole exchange within the framework of envelope 
approximation. We note that Takagahara's condition [16|] for the vanishing of LRE in 
spherically symmetric quantum dots constitutes only a sufficient condition for quenching of 
LRE. It does not explain, for example, why LRE vanishes in disk-like QDs with rotational 
(Coon) symmetry. These QDs are "squeezed" in one direction and the ground state envelopes 
of single-particle states will contain angular momentum components that are higher than 
Yqo(8, <p). Motivated by recent experiments, we will investigate, within our model, the effects 



3 



of the external electric field and size-scaling of LRE interaction. 

We derive here effective "four band" excitonic Hamiltonian which includes effects of LRE 
interaction. The elements of the effective Hamiltonian are expressed explicitly in terms of 
envelope functions. The "microscopic parts" of single-particle orbitals are integrated out and 
enter the effective Hamiltonian as numerical parameters. The number of bands is truncated 
to two conduction and two valence bands for the simplicity of the interpretation. We present 
explicit expression for the matrix element responsible for the "bright" exciton splitting and, 
therefore, establish a new sufficient condition for LRE quenching. 

An excitonic fine structure for a model system in which the confining potential has a 
form of two-dimensional-like anisotropic harmonic oscillator (2DLAHO) is considered as a 
function of lateral anisotropy. We present an explicit expression for the "bright" splitting 
of excitonic ground state as a function of lateral anisotropy. It is found, in agreement 
with previous work [201 ]. that, within the envelope approximation, the "bright" ground state 
exciton splitting vanishes in the case of laterally isotropic confining potential. The quenching 
of "bright" exciton splitting coincides with vanishing of the matrix element responsible for 
the "bright" exciton splitting in our effective Hamiltonian. 

We also analyze the effect of the lateral electric field on the excitonic fine structure of 
our model quantum dot. We find that the "bright" exciton splitting decreases due to the 
spatial separation of electron and hole envelopes. 

Finally, the scaling of the "bright" exciton splittings with system's size is analyzed. It is 
found that the scaling of "bright" exciton splittings differ from the laws established using 
simple dimensionality arguments. 

Effective units of length and energy are used throughout unless otherwise specified. The 
lengths are measured in effective Bohrs a e B = eh 2 /(m*e 2 ), where e is dielectric constant and 
m* is the conduction band effective mass. Energies are measured in effective Hartrees, 1 H = 
m*e 4 /(eh) 2 = 2 Ry*. For example, using material parameters for GaAs, a e B = 97.9 A and 
I H = 11.86 meV. 
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II. THEORY OF EXCITON FINE STRUCTURE IN ENVELOPE FUNCTION AP- 
PROXIMATION 



We now describe the single particle states and the exciton fine structure in the envelope 
function approximation. The single-particle orbitals of the electron in a quantum dot are 
two-component spinors (two-row columns) written as 






, 1/5) = 1 


: 









a (r)|a> + 6 (r)|/?>, |a> = ,!/?>= • (!) 



A dagger sign t will denote complex conjugation for complex quantities. 

With and (c and h) electron and hole creation (annihilation) operators, the inter- 
acting electron-hole Hamiltonian can be written as 

H x = H e + H h + H int , H e = Yl £C i c U, H h = J2 ^h]hj, 

» 3 

H int = J2 C i kt Ml (~V&l + Viklj) > 
ijkl ^ 6 ' 

v, m = V(i,r,k,i) = / /^ ^W^") (2) 

j j |ri - r 2 | 

where ff e , and are electron, hole, and electron-hole interaction Hamiltonian, re- 
spectively. The electron-hole interaction Hamiltonian consists of two parts: V c /e is the 
direct electron-hole Coulomb attraction screened by static dielectric constant e, V E is the 
electron-hole exchange interaction. The question of screening of electron-hole exchange is a 
subtle one and it is generally agreed that, at least, parts of the electron-hole exchange should 
be screened. In Ref. 27| (page 252), the long-range electron-hole exchange interaction in 
reciprocal space is screened by the high-frequency dielectric constant. We will assume, for 
now, that our exchange interaction V E contains screening implicitly. 

To obtain excitonic states, Hamiltonian (j2j) is diagonalized in the basis of all electron- 
hole pairs of the type vhJ\g.s.), where \g.s.) denotes a many-body state with fully occupied 
valence and empty conduction bands. 

In this work, the hole and electron single-particle states are computed in the effective mass 



approximation (EMA) 



27 



28j | which neglects band coupling in the single-particle states. In 



EMA, the hole 4>h and electron <p e single-particle states are uniquely specified by band label 
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(valence vj' z and conduction ca') and envelop index (r and p), for example, 



M r ) = F vj' z ( r )uvj' z (r), 
e (r) = ^(r)vW- 

(3) 

Here, u(r) is the periodic part of Bloch eigenstate at k = 0. In what follows, u(r) is a 
two- component spinor. We will assume that "bra" electron-hole pair is described by indices 
cas and vj z q, respectively. 

The excitonic Hamiltonian matrix element in the basis of electron-hole pairs is 

(g -S.\h v j zq C cas H X c \ a i p h\j^ T \g -S.) = $ v j : ,q ) vj z rfic(TS,C(T'p{.£ v j z q + e C as) 

-V c (cas, vj' z r; vj z q, ca'p)/e + V E (cas, vj' z r; ca'p, vj z q) (4) 
The effective Hamiltonian which involves only envelope functions is obtained by inte- 



grating out "microscopic" degrees of freedom (u). The derivation follows that of Ref. [20 ] . 
Assuming symmetry of the crystal, two conduction band "microscopic" functions \cl/2) 
and |c — 1/2) with spin projections s z = 1/2 and s z = —1/2 are written as 

Mcl/2 (r) = (r|cl/2) = xs(r) I ^ I , w c _ 1/2 (r) = (r|c - 1/2) = Xs (r) I ° I . (5) 

In Eq. ([5]), we assumed that "microscopic" conduction band function are of pure s symmetry. 

Two valence band "microscopic" functions \v — 3/2) and \v + 3/2) with hole angular 
momentum projections j z = —3/2 and j z = 3/2 are taken as 

i Y n(r)\ ( 

WB-3/2(r) = (t\v - 3/2) = x P (r) , u v+3/2 (r) = (r\v + 3/2) = x P (r) 

/ lii-i(r) 



(6) 



In Eq. we assumed that "microscopic" valence band functions are of pure p symmetry. 
The eigenf unctions of angular momentum Xp( r )^ii(^) an d Xp( r )^i-i(^) can De expressed in 
terms of Cartesian functions p x and p y . We, therefore, are neglecting p z contribution to the 
"microscopic" valence band functions. 

In what follows, we give a brief derivation of the effective excitonic Hamiltonian. 
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A. Calculation of Coulomb Direct Matrix Elements. 



The unscreened Coulomb direct matrix element is given by 

T/C7 " ■ M f ( a a <7i(ri)g 2 (r 2 ) 

V (cas,vj z r;vj z q,ca p) = / / dridr 2 - 



ri - r 2 



q 1 (r 1 ) = F^( ri )«L(ri)^(ri)^(rO 
g 2 (r 2 ) = Fj,(r 2 )^,(r 2 )F4(r 2 )^(r 2 ) 



(7) 



where we explicitly presented electron and hole single-particle orbitals as a product of an 
envelope F and a "microscopic" part u. The unscreened Coulomb attraction matrix ele- 
ment is just Coulomb interaction between two "transition densities", where each "transition 
density" is a product of two electron orbitals or two valence orbitals. The matrix element is 
approximated as 

V c (cas, vj'r; vj z q, ca'p) = <5 W <W / / dr^ ^ 1 ^* 1 ^* 2 ^^ . (8) 

J J I r l — r 2 1 

B. Calculation of Coulomb Exchange Matrix Elements. 

The exchange matrix element can be written as 

V (cas,vj z r;ca p,vj z q) = / / dridr 2 - 



ri - r 2 



q 1 (v 1 ) = F^ul^F^uM, 
g 2 (r 2 ) = F c p cr ,(r 2 )M C(T /(r 2 )F^ 1 ,(r 2 )M|; j ,(r 2 ). 

(9) 

The exchange matrix element can be thought of as a Coulomb interaction of two "transi- 
tion densities" , where each "transition density" is a product of electron-hole single-particle 
orbitals. 

We decompose exchange integral (jHD into the short-range and long-range contributions 
in real space. We will refer to the whole integration region as Born-von Karmen cell (BvK 
cell) and to the individual unit cell within BvK cell as Wigner-Seitz cell (WS cell). The 
double integration over BvK cell (which consists of N ce u WS cells) is replaced by N cea x N ce u 
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integrals over WS cells 

N 11 N 7i 

p p A *cell p p ly cetl p 

/ / dridr 2 — ► / / dri<ir2+ V] / 

8, 7 = 1 



dr 1 dr 2 , 



,3 = 
i ^ 3 



GWS(Ri) Jr 2 eWS(R.j) 



(10) 

where Rj, R., label positions of the WS cells. The first term in (|TU|) consists of N ce u integrals 
in which and r 2 go over the same cell, whereas the second term consists of N ce u x (N ce u — 1) 
integrals in which ri and r 2 go over two distinct WS cells. The first sum in which ri and r 2 
go over the same cell is the short-range exchange V^ R , whereas the second sum in which r± 
and r 2 go over two different WS cells is the long-range exchange V^ R . 

After some algebra which involves multipole expansion of l/|ri — y 2 \, we obtain the 
following expression for the exchange matrix element 

V E (cas, vj' z r; ca'p, vj z q) = Vg R (cas, vj' z r; ca'p, vj z q) + V^ R (cas, vj z r; ca'p, vj z q), 
V s E R {cas,vj' z r;ca'p,v Jz q) = E SR {Hf^, J drF^r)F^(r)F^(r)F^(r), 



V E R (cas,vj' z r; ca'p,vj z q) = ~y^ 2 (( d ! 



• Y 

cavjz ) 



d ^) J dvF^(r)F v l(v)F^(v)F:l(r) 



3 

^ 2 E( d c^j;(d^) 5 / / 

7,5=1 J J 



drjdrl 



ri - r2 



-dr 1 dr 2 



(11) 



where E$r and /i 2 are two numerical constants parameterizing short- and long-range ex- 
change interaction, respectively, (Hgg)^^, is an element of short-range exchange "micro- 
scopic" matrix 



H 



int 
SR 



/ 
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(12) 
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and d° C(TV j z are "microscopic" dipole elements 



fd° ■ ) (d° ■ ) (d° ■ ) 

V^-arvjzlx V^-ccrvjzJy V^cavjzJz 



\c(?,vj z ) ||,-|) 


-1 


— i 







\ca,vj z ) ||,+|) 











(13) 


\c(T,vj z ) | - §,-§) 













\ccr,vj z ) | - §,+§) 


1 


—i 


0. 





The numerical parameters E$r and fi 2 can be determined as to reproduce the excitonic 
fine structure in bulk semiconductors and, therefore, implicitly contain screening effects. 

Taking into account ( fl2i) and i fTBl . a "block" of Hamiltonian corresponding to the ex- 
change interaction between two electron-hole pairs F£F£ and F*F£ can be presented in the 
form 

/ 









|ax',^) 


\ca\ vf z ) 


\ca',vj' z ) 






II _3\ 
1 2' 2/ 


II +3\ 
12' ~ll 


1 - I -3\ 

1 2' 21 
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1 2' '2/ 


\ccr,vj z ) 


II _3\ 
1 2> 2' 
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d 21 








X SRE,L . X LRE,L . X LRE,N 

d +d +d ) 



(14) 

where we separated different contributions to the exchange based on their origin (SRE or 
LRE) and integral type ("local" and "nonlocal"). The contributions are 



S S RE ' L = E SR j dvF?(r)F*(r)Fv(r)F?(v), 



fiLRE,N _ _^2 ^> 



■yy) 

5 lre,n = v^(l/2 S) 3/2r; -l/2p, -3/2q) = fx 2 (R xx - 2iR xy - R 
6 lre,n = v^(_l/2s, -3/2r; l/2p, 3/2q) = li 2 (R xx + 2iR xy - R 



yy) ■> 



yy) > 



(15) 



With regard to the above exchange expressions we note the following: (a) The short- 
range exchange causes splitting between "bright" (|^,— |), | — |,+|)) and "dark" doublets 
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(||, +|), I — §, — §}) by moving "bright" doublet up in energy by 8q RE,l . SRE does not split 
the "bright" doublet. The integral describing SRE is "local" in nature and involves overlap 
between two electron-hole pairs, (b) Long-range exchange contains expressions of two types: 
a "local" term 8q RE ' l which arises from fi 2 ((d° rt)(7 )"l" ■ d° rW ) and "nonlocal" terms 8 LRE,N , 
8 12 RE ' , and 821" ' which involve differentiation operators applied to electron-hole pair 
envelops. The "nonlocal" terms describe dipole-dipole interaction between electron-hole 
transition densities, (c) The splitting between "bright"- "dark" states is affected by LRE 
through the "local" LRE term 8 LRE ' L and "nonlocal" LRE term 8 LRE,N . LRE, in general, 
splits the "bright" doublet by coupling electron-hole pairs with "completely opposite" z- 
projections of angular momentum, for example, | — |, | ) and |) • The terms which 
are responsible for the "bright" exciton splitting are 8 E 2 RE ' N and 82 RE,N - the "nonlocal" 
expressions arising from long-range exchange interaction, (d) The splitting within the "dark" 
doublet is not described in our model. 

The splitting of "bright" exciton levels vanishes provided that 8^ RE ' N and 82 RE ' N vanish 
which constitutes a condition for LRE quenching within our model. In the next section, 
we demonstrate numerically and analytically that this condition is fulfilled in the case of 
isotropic 2D HO confining potential. 



III. EXCITON FINE STRUCTURE OF QUASI- TWO-DIMENSIONAL 
ANISOTROPIC HARMONIC OSCILLATOR (2DLAHO) QUANTUM DOT 

In this section, we will investigate excitonic fine-structure for a model quantum dot in 
which the confining potential is of 2D-like anisotropic harmonic oscillator type. We note 
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that HO oscillator spectrum has been observed in self-assembled quantum dots 
The EMA equations for holes are 

Hhh = ~7yj\7 ( V * + V ?) + \ M W^ Khh* 2 + 4mV 2 ) - ^TF^ V ' + \ M ^\ hh z\ 
*M\\,hh A ZM±h h Z 



U x ,hh — ^hhi^ + £)) ^yM — hH , ^x,hh U yM — {^hh) 2 ~ COUSt, U Z: hh ^hhi 



l+t 

(16) 



where M\^ hh and M^ hh are effective masses, cj 7 (7 = x,y,z) denotes confinement frequency 
in x, y, and z directions, respectively, determines the typical energy scale of the confin- 
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ing potential and, together with the mass, the confinement. The energies and lengths are 
expressed in the effective units. 

EMA equations for electrons have a similar form 

H ee = ~\ {V 2 X + V 2 + V 2 ) + ~ « ee x 2 + ul ee y 2 + uj\ ee z 2 ) , 



w2 



^x,ee = ^eei 1 + *)> U V,ee = 7~7> ^x,ee^y,ee = {^leY = COHSt, L0 z , ee > 

-L ~ r C 



(17) 



The excitonic fine structure is studied as a function of anisotropy t. For example, for 
t = —0.5, uj x = (l/2)u° and u y = 2a; , the confinement along x is weaker than along y. For 
t = 0.5, ct^ = (3/2)u; and u; y = 2/3a; , confinement along x is greater than along y. In the 
case t = 0, the confining potential for holes and electrons is isotropic in the lateral direction. 
By changing anisotropy t, we control the shape of the confining potentials for holes (fTBl) and 
electrons ffTTI) and, as a consequence, the single-particle energy spectra and eigenfunctions. 

We have assumed that the confinement in the vertical direction z is much stronger than 
confinement in the lateral direction uo z ^> u°. Therefore, we will always assume the ground 
state solution in z direction and the single-particle energy spectra for holes ffTBl) and electrons 
(EE7D is 2D-like 

+ 2~) Ux+ ( m+ 2~) Uy + \ Uzj ^ 
where n and m 2DLAHO quantum numbers. In what follows, we use numerical parameters 

M\\ thh = 9.930853, M ±M = 5.218662, M Ue = 1.0, M ± , ee = 1.0, oj z>hh = 20.155737, and 

u z ^ ee = 105.185977. 

The eigenfunctions are products of one-dimensional Hermite polynomials and exponential 
function 

^nmo(r) = 4> n (x)4>m(y)4>o(z), (19) 

where the eigenfunction in one of the Cartesian directions (for example, x) is given by 



Wn{.X) = Y ^7 \~^~) eXP \ 2~ X ) H nW muJ xX) 



^,(x) = (-irexp(x 2 )^ r exp(-x 2 ) 



(20) 



Figures (la), (lb) and (lc), (Id) show single-particle spectra for electrons and holes, respec- 
tively, for two different characteristic energy scales uj° . Solid horizontal line (around 52.6 
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and 10.1 H for electrons and holes, respectively) corresponds to the confinement energy in 
vertical direction uj z /2. 

The single-particle energy spectra are given by (ITS)) . The single-particle ground state is 
the nodeless s envelope which corresponds to n = m = 0. The "excited" levels are organized 
in shells referred to as p, d, f, and so on with characteristic "spatial" degeneracies of 2, 
3, 4, etc. in the case of isotropic (t = 0) confinement potential. The dispersion E(t) of 
s energy level with anisotropy is weak compared to the dispersion of "excited" envelops. 
Strong anisotropy might lead to the level crossing where the single-particle energy of <i-type 
envelope is below the single particle energy of p-type envelop. The level crossing happens, 
for example, at t pa —0.3 for the second excited state on Fig. (la). Due to the smaller 
effective mass of electrons, the spacing between electron single-particle levels is larger than 
that of the hole levels. 

Figures (le) and (If) show "noninteracting" electron-hole and "interacting" exciton en- 
ergies as a function of lateral anisotropy t for two different characteristic energy scales u°. 
The black squares on Figures (le) and (If) correspond to the "noninteracting" electron- hole 
pair energies S vaqtVa ' r 5 CTStCT ' p (e^ acl + e e CTS ). The empty circles are obtained by diagonalizing 
the full excitonic Hamiltonian. Each of the empty circles is actually a multiplet of 4 excitons 
with fine structure determined by the electron-hole exchange interaction. The "noninter- 
acting" electron-hole and excitonic spectra look quite similar. The Coulomb electron-hole 
attraction simply decreases the exciton energy. The Coulomb attraction-induced mixing of 
electron-hole pairs in an exciton is small due to the large separation of single-particle levels 
on the scale of magnitude of screened Coulomb attraction. 

In the lower part of the "interacting" excitonic energy spectra and for small anisotropies 
\t\ < 0.2, the dispersion Ex{t) strongly resembles that of the dispersion of single-particle 
levels. This happens because Lo\ h < u>® e and the lowest energy electron-hole pairs follow the 
order of hole levels s e Sh, s e ph, and s e dh- 
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A. LRE interaction for s-type envelopes 

Consider the case when exciton is given by the product of ground state envelopes of 
electron and hole single-particle states (s e Sh exciton). In this case, 

F^(r)F v (r) = N c N v exp (-a x x 2 - a y y 2 - a z z 2 ), 

N = (2a^ 1/4 ^2ay^ 1/4 (2a^ 1/4 ^ = ^2a xv ^ 1/4 ^2a yv ^ 1/4 ^2a„^ 1/4 

Qxc ~t~ (%XV7 Q-y &yc ~t~ Qyvi &z &zc ~t~ &zvi 
_ ^x,ee _ L*Jy,ee _ ^z,ee 

(%xv 2 ' CXyv 2 ' ^ zv 2 ' 

(21) 

where N c and N v are normalization constants, a lc and a lv (7 = x, y, z) denote confinements 
of electrons and holes in x, y, and z direction, respectively, and a 7 = a lc + a lv denote 
confinements of electron-hole envelope. 

The matrix elements responsible for "bright" exciton splitting for s e Sh exciton are given 

by 

^12 = A* (-R-xx 2iR xy Ryy) ) 
^21 ' = A 4 {R-xx "I" 2iR xy Ryy) -i 

(22) 

where 



J J \ ax{ J |ri — r 2 | a x a y a z 

J ° J ° (^^(icos^+^sin 2 0) + ^cos^) 

(r 2 )^ f (r 2 ) ^ (Ay^) 2 ^ 

: : — aridr 2 = n^irl v {a x , a v , a z ), 

|ri-r 2 | a x a y a z 

f w/2 f w/2 sin 3 6 sin 2 <f)d6d<f) 

l y {a x ,a y ,a z ) — / / — - - \3/2' 

J0 J0 (sin 2 * (i cos^ + sin 2 0) + £ cos^ 6>) 



4(« 
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It follows from fl23|) that S 12 ' and <5 21 ' vanish provided that the "confinement" of 
electron- hole pair envelope is identical in x and y directions (a x = a y ). Figure 2 illustrates 
this. Figure 2 shows "bright" exciton doublet splitting as a function of lateral anisotropy 
t. The doublet splitting energy is defined as E y x — E x , where E x and E v x are ground 
state "bright" exciton energy levels polarized along x and y directions, respectively. For 
t < {uj x < uj y ), E y x > E x and the "bright" ground state is dipole active along the x- 
axis, whereas for t > (uj x < uj y ) the "bright" ground state is dipole active along the 
y-axis. In our model, the "bright" ground state is always dipole-active along the axis of 
weaker confinement. The insert to Figure 2 shows R xx and R yy "nonlocal" contributions to 
LRE exchange as a function of anisotropy computed from equations (12"51) . We can see that 
R xx increases and R yy decreases in magnitude as t goes from -0.5 to 0.5. The insert also 
shows the difference R xx — R yy which determines "bright" exciton splitting as a function 
of anisotropy. We can see that R xx — R yy steadily decreases as a function of t and passes 
through zero at t — 0. t = corresponds to the laterally isotropic confining potential. In 
this case, R yy = R xx and the exchange matrix element which couple two "bright" exciton 
levels vanishes. As a results, the splitting between "bright" excitons vanishes as well. 



B. Application of lateral electric field 

In the previous section, we have demonstrated that the magnitude of the "bright" exciton 
exchange splitting can be controlled through the shape of the confining potential. In practice, 
one may, for example, try to quench the "bright" exciton splitting by picking "symmetric" 
dots from a large number of samples grown under different conditions [fj. Nethertheless, 
the growth of QD is essentially a random process and control through the QD shape is hard 
to achieve. 

It is, therefore, of great interest to examine the effects of external fields on the excitonic 
fine structure. For examp.e, QDs cau be p.aced Q betweeu Schottky gates for the application 
of vertical and lateral electric fields. Recently, it has been demonstrated theoretically [30] 
and experimentally [12j that by applying an in-plane electric field it is possible to fine-tune 
photon cascades originating from recombination of multiexciton complexes in QDs. 



The electron-hole exchange interaction was treated in Ref. 30J using empirical exchange 



Hamiltonian. This subsection discusses the effects of the lateral electric field on the exchange 
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matrix elements of our effective EMA Hamiltonian (fl4j) . We will focus on the "bright" 
exciton splitting as a function of the lateral electric field. 

Within our model, application of lateral electric field F = Fe x + Fe y of magnitude \F\ 
displaces the "origin" of electron and hole envelops in xy-p\a.ne from (0, 0) to (x^, y$) and 
(xq,t/q), respectively. The single-particle energy levels are rigidly shifted by the Stark shift 
whereas the spacing between the single-particle energy levels of HO is not affected. The 
"separation" of electron and hole envelopes in xy-plane is determined by equation 

Ax = x e -x% = eF\ + 1 2 ) , 

Ayo = y e -yZ = eF[-^+ 1 ) . 

\^y,ee M \\-hh^hh J 

(24) 

Evaluating the integrals which control "bright" exciton splitting R XX (F), R yy (F), R xy (F) 
as a function of field F for s e Sh electron-hole envelope, we obtain 

RUF) = R xx (0) exp (-^^Axl - ^yiAy^) , 

R yy (F) = R yy (0) exp f_ 2 _^l Ax l _ ^fviAyl) , 

R xy (F) = (25) 

where R xx (0) and R yy {0) are integrals (123!) in the absence of the electric field (F = 0). Field 
dependence of R XX (F) — R yy (F) which determines the "bright" exciton splitting is shown on 
Figure 3 for five different initial values of lateral anisotropy. At F = 0, the "bright" exciton 
splitting attains maximum value which is determined by the initial shape anisotropy of the 
confining potential. The larger the magnitude of the anisotropy, the larger is the initial 
(F = 0) "bright" exciton splitting. As F increases in magnitude, the electron and hole 
envelops are pulled out in opposite direction and the magnitude of the splitting is reduced. 
Since the separation between the envelopes depends on F 2 , the splitting is independent of 
the sign of the field F . 

It is interesting to note that the "rate of quenching" of the splitting depends on the 
initial anisotropy - the stronger the confinement, the larger is the drop in the magnitude 
of the splitting. One can use lateral electric field to produce two identical "bright" exciton 
splittings for two dots with different initial anisotropies. This happens, for example, at 
|F| « 0.6 for two dots with initial anisotropies of t = —0.5 and t = —0.4, respectively. 
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C. Scaling of "bright" exciton splitting 



The question of size-scaling of exchange interactions in nanosystems has received a lot of 
attention, particularly, in connection with the so-called Stokes ( "red" ) shift of resonant PL 



spectra with respect to the absorption edge (see, for example, Ref. (18]). In this section, we 
examine the size scaling of "bright" exciton splitting as a function of system's size. 

Suppose that R is a characteristic size of a 2D-like system. From the normalization 
condition 

rfr|F c (r)| 2 = l, I dv\F v {v)\ 2 = 1, (26) 



envelope functions scale as 1/R 2 . Based on this "normalization" argument, the nonlocal 
LRE integrals R xx ,R yy , and R xy that determine "bright" exciton splitting scale as 

wwhw R2R2 K w < 27 > 

The dependence ( 1271) is expected to hold in the strong confinement regime. Of course, one 
has to keep in mind that the exchange matrix element responsible for the "bright" exciton 
splitting is a linear combination of nonlocal LRE integrals R xx , R yy , and R xy . Therefore, 
the actual size dependence of "bright" exciton splitting might be different from 1/R 3 and 
depend strongly on the "nature" of the envelop functions involved in the LRE integrals. 

Suppose the HO potential is of the type Vho{. x iD) — V (x 2 + y 2 )/R 2 , where V is some 
characteristic energy scale and R is characteristic size determined by HO frequency uj°, mass 
Mm, and characteristic energy scale Vq 



R =hpk= e i L < 1 =Jt£p L =\jm- (28) 

where I and L are "confinement" lengths. If Vq is kept constant, R cx l/(M||Cc; ). We will, 
therefore, calculate the "bright" exciton splitting as a function of l/(M||u; ) keeping the 
ratio of electron and hole confinement frequencies constant uj^/uj^ = const and the lateral 
anisotropy t fixed. Since ^ee/^hh = cons t by increasing/decreasing uP ee we automatically 
increasing/decreasing u)\ h . 

Figure 4 shows the scaling of the splittings for the ground and excited "bright" exciton 
levels with size. The lateral anisotropy is fixed to t — —0.1 in all the calculations. We 
find that splittings decay as the size of the system increases. We find that the ground state 
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exciton splitting is size-insensitive whereas the splittings of the excited bright excitons scale 
oc l/R 13 . This scaling is different from l/R? dependence expected from the "normalization" 
arguments ([27]) . The deviation might be due to the fact that the exchange matrix elements 
that determine the "bright" exciton splitting 5 12 and 621 involve a linear combination of 
nonlocal exchange integrals R xx , R yy , and R xy . Moreover, R xx and R yy enter the expression 
for 5\2 and 821 with different sign and some cancellation of terms may occur. Therefore, 
the size dependence of "bright" exciton splittings may depend strongly on the details of the 
excitonic wavefunction. 

IV. CONCLUSIONS 

A four-band electron-hole excitonic Hamiltonian is derived within the effective mass ap- 
proximation (EMA) which takes into account the electron-hole exchange interaction. The 
matrix elements of the effective Hamiltonian are expressed explicitly in terms of electron and 
hole envelopes. The "microscopic" parts of single-particle orbitals are integrated out and 
enter our Hamiltonian implicitly through numerical parameters. The matrix element respon- 
sible for the "bright" exciton splitting is identified and analyzed. An explicit expression for 
this matrix element in terms of electron and hole envelopes is presented. An excitonic fine 
structure for a model system with 2D-like anisotropic HO confining potential is considered. 
It is found that in the case of isotropic potential, the "bright" exciton splitting vanishes. 
Within the formalism of our effective excitonic Hamiltonian, the effects of the lateral electric 
field on the excitonic fine structure were considered. It is found that the excitonic structure 
can be tuned with lateral electric field and that the magnitude of the exchange splitting is 
reduced but never to zero. The origin of this reduction is the separation of electron-hole 
pairs. Finally, size- dependence of the "bright" exciton splittings for the ground and excited 
state excitons was investigated and it was found that the size scaling of "bright" exciton 
splitting can be different from the laws established using normalization conditions. 
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FIG. 1: Single-particle and Exciton Energy Levels as a function of lateral anisotropy t of the 
confining potential. 
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FIG. 2: "Bright" exciton doublet splitting as a function of lateral anisotropy t of the confining 
potential. Insert: R xx (squares), R yy (circles) contributions to "nonlocal" LRE matrix elements as 
a function of anisotropy of the confining potential; R xx — R yy (stars) difference of two "nonlocal" 
LRE contributions that determines the "bright" exciton splitting. For the isotropic confining 
potential t = and R xx — R yy = 0. 
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FIG. 3: "Nonlocal" LRE R xx — R yy which determines "bright" exciton splitting as function of 
lateral electric field F. At F = 0, R xx — R yy attains maximum determined by the initial anisotropy 
of the confining potential. As magnitude of field F increases, R xx — R yy decreases due to the 
separation of electron and hole envelopes. 
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FIG. 4: Size scaling of the "bright" exciton doublet splittings. Lateral anisotropy is kept constant 
t = —0.1. The doublet splittings are plotted as a function of l/(Mn ee uj® e ). The ratio uj® e /LU® h is 
kept constant. The ground exciton is of s e — Sh type. Two excited "bright" exitons correspond to 
s e — dh electron-hole pairs. The solid lines are fit to the power law C/R n , where C is a constant, 
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